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t> ! INTRODUCTION 

^3 | Recently there has been a lot of papers on matrix models and superstrings, induced 

by the work of Banks, Fischler, Shenker, and Susskind (hep-th/ 9706168). I refer to 
■ Makeenko's talk at this meeting for a general review of this subject. 

Most of the workf] reported in this talk has been done together with Fayyazuddin, 
Makeenko, Smith, and Zarembo 1 . As explained in Makeenko's (virtual) talk at this 
meeting, we started from the work by Ishibashi, Kawai, Kitazawa, and Tsuchiya 2 , who 
proposed that type IIB superstrings in 10 dimensions are described by the reduced 
action, 

Sikkt = a (~Tr [Ap, A v ] 2 - ^Tr {^[A^]) 

where A^ and ipa are n x n matrices. A sum over n is implied, with weight exp(— (3n). 
Later the sum over n has been replaced by a double scaling limit 3 . 

In our paper 1 we discussed various problems associated with eq. (1), and we 
proposed a different model with action 

Snbi = -^Tr (y- 1 ^, A,] 2 ) + V(Y) - |Tr (^[A,^]), (2) 

where the potential is given by 

V(Y) = (3 TrY + 7 Tr In Y. (3) 



*The paper in reference 1 has been published in Nuclear Physics B. Unfortunately, the editors of that 
journal used an early draft of the manuscript, which contained several typhos. For this reason I 
cannot recommend the published version, but refer to the version in the Archives. 
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The partition function is thus given by 



dA^dipdY e~ 5NBI . (4) 

We select the constant 7 in such a way that the result of the Y-integration is as close 
to the superstring as is possible. This turns out to mean 

1 

l = n- -, (5) 

as we shall see later. 

Physically the model S^bi is motivated by a GUT scenario: Suppose one has a 
field theory valid down to the GUT scale. Then, in our model the group is SU(n), 
with n large. As we shall see, this type of GUT model then leads to superstrings if 
n — > 00. For n finite, supersymmetry is broken, as is expected for energies below the 
GUT energy. Thus, superstrings can emerge from a GUT type of model. Of course, 
the model with action ,S NBI is not a realistic GUT model. 

Under the SUSY transformations 

5*P = l -{Y-\ [Ap, A u ]}T^e, 5A, = iel>, (6) 
the action transforms like 

5S NBl oc e ^i...A 7 Tr ^ m (Y°Y u T Xl ...T X7 ) mp e p {[A a , Ap], [A^Y' 1 }}) . (7) 
It can be shown that 

£5nbi — ► for n — > 00, (8) 

so the action is supersymmetric in the limit n — > 00, but for finite n the symmetry is 
broken. 



THE Y-INTEGRATION 

The integration over Y can be done exactly. Consider 

F{z) = JdY exp (—^Ty(Y~ 1 z 2 ) — (3 TrY — 7 Trln Y^ , z 2 = -[A^,A U ] 2 . (9) 
The "angular" integration is of the Itzykson-Zuber type, so we get 

Hz) = const. J] J ^ A(1/y) A (z 2) e"«&^^ h * (10) 

Here the z^'s and y^s are the eigenvalues, and 

A(ar) = Y[(xi - X j) = det x k ~ l (11) 

i>j 

is the Vandermonde determinant. We only integrate over the positive eigenvalues of Y. 
Thus we get 

/•oc 

A(z 2 )F(z) = const. / Y[dyi y?- 1 JJ(^ - y 5 ) e ' a ^' z '^-^' y ^^ hlv ' . (12) 
J ° i l>3 



2 



This can be rewritten as a determinant 



A(z 2 )F(z) = const, det / ~ky k ~ l e 



dy __ k _i az 1/4y-(3y 



const, det 

ki 



. f) k - 1 / PtF „ m 

(13) 



i)k-i dk 1 / r e -v : >■-. 



Here Zj is by definition the positive square root of zf. This determinant can be evaluated 
using basic properties of determinants, and the result is 1 

A(z 2 )F(z) = const. A{z) (14) 

This result is exact, and hence it is valid for any n. 

The sum over eigenvalues in the exponent has the following interpretation, 

where we used an integral representation of the square root. Thus, 

I> = ^ W Tr exp(-tK,^] 2 ) = Tr ^j-\A^A v \\ (16) 

valid in Euclidean space. 

The partition function therefore becomes 

Z = J dA^dipdY exp (^Tr(Y^[A^ A u f) — (3 TrF — 7 Tr toy) 

= const. J dA £* exp (-y/tf Tr^/-[A„ Arf - |Tr (#^,^])) . 

(17) 

This is the exact result of the y-integration. In order to get the square root it is 
important to use the value of 7 given in eq. (5). 

This result can be expressed in an alternative form, at the cost of introducing an 
auxiliary Hermitean field M. We use the identity 



U^M + Zj) 



t<3 

to obtain 

1/2 

12 ' 



const. Vdetz~ J dMe~ Tr zM \ (18) 



Z = const. J dA^dipdM (det yj-[A^A v ] 

xexp (-Tr ((y/tf + M 2 ) ^ -[A„ Atfj - |lV {^[A^]) ) . (19) 



The field M is essentially trivial, with a "classical equation of motion" M — 0. 
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ON THE WEYL REPRESENTATION AND THE APPROACH OF THE 
COMMUTATOR TO THE POISSON BRACKET 

The square root occuring in the result above is somewhat reminisent of the Nambu- 
Goto square root. If we could replace the commutator in the square root in eq. (17) by 
the corresponding Poisson bracket, we would have a partition function which is very 
similar to the one for the superstring. 

This problem has been discussed by Hoppe 4 , and in different settings by a number 
of other authors 5 ' 6, 7 . It turns out that making some assumptions, one has the limit 

" [A,B] -+i{A,B} PB ", for n ^ oo. (20) 

Here {, } PB denotes the usual Poisson bracket. 

We refer to the literature for a detailed discussion. Here we shall follow Bars 7 , and 
consider a torus (although this restriction is probably not important 8 ). The case of a 
sphere was discussed in ref. 5. A Hermitean matrix can be expanded in a Weyl basis, 

(4.)5 = CiX>i('k)J, with k=(k u k 2 ). (21) 

k 

The matrix Z k can be expressed in terms of the n x n (n =odd) Weyl matrices h and 
g, which satisfy 

h n = g n = 1 and gh = ujhg, with u = e 4ni/n . (22) 
The explicit form of these matrices are 

h = diag(l, u, u 2 , ...,u n ~ v ) 

<7} = <Jj +1 , * = 1,2,.. g? = except for g n x = 1. (23) 
The SU(n)-generators Z k are then constructed as 

Z k = — c/ lfc2/2 /iV 2 , (24) 
Att 

since the powers of h and g are linearly independent for ki,k 2 = l,2,...,n — 1, are 
unitary, close under multiplication, and are traceless. Using that 

Tr />V 2 = n5 ku0 5 k2 , , (25) 

we easily see that 

Tr l k l p = j^yS k+p , . (26) 
Thus the expansion coefficients in eq. (21) are given by 

% P = KY = §^Tr (l p A,). (27) 



Also, from the relation 

we get 



g k lh k 2 =UJ k 1 k 2h k 2g k 1 ^ 



[Ip, k] = sin ^^-p x kj Z p+k , 



(29) 



where 

p x k = - hP2- (30) 
Using the expansion (21) we get 

[4» aja = c? E < a "^ sin ( v p x q ) (/p+q) * • (31) 

p,q ^ ' 

This can be compared with the similar expression for the string variables X M (cr, r), 
where we have the expansion 

X^a, t) = C 2 E a ™ exp( y iam 1 + irm 2 ), (32) 

m 

leading to 

{X„ XJpb = C\ « (P x q) exp(*(p + q)<r). (33) 
p,q 

Now it is clear that if 



lim V ... = V lim (34) 

n—>oo ' J ' J n^oo 

modes modes 

then we have by use of (26) 

Tr [Ap, A u ] 2 -> const. J dadr{X^X u } 2 for n -> 00. (35) 

It is obvious that the commutativity (34) is only valid if the infinite modes are unim- 
portant. This is, however, not true e.g. for the bosonic string. If we fix the end points 
of this string at some distance, then there is a critical distance (essentially the inverse 
tachyon mass) at which the string oscillates so wildly that this behaviour can only be 
reproduced with an infinite number of modes. Below this distance the "string" becomes 
a branched polymer, and hence is no longer a string. 

For superstrings this problem does not arize, and hence there is at least no obvious 
reason why the limits cannot be intechanged as in eq. (35). In the following we assume 
that eq. (35) is correct for type IIB superstrings. 

Since we are interested in the square root of the squared commutator, the result 
(35) is not enough. Using eq. (26) and repeated applications of the formula 

IrJr = (n/An) exp(27ri(m x r)/n))Z m+r (36) 

one can easily derive 

Tr / mi / m2 .../ ms -> n(n/4ir) s 5 mi +m 2 +...+m s ,o, for n -> 00 (37) 
to leading order in n. Using this in eq. (16) we obtain 

I 1 r {0+) dt 

Tr yJ-lA^Atf = -£j= J ^ ^Tr exp(-t[A„, A u f ) 

1 [ {0+) dt f 
~^ 4~7^ J tsTs J d(Tdrex P(t const. {X^,X u }l B ) 

= const. J dadr^J {X^, X U } PB for n — > 00. (38) 
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Thus we see that in the leading order the Nambu-Goto square root arises as the limit of 
the square root of the corresponding commutator. However, it should be remembered 
what was said before about strings with tachyons. They do not allow the interchange 
of limits as in (34), and hence the result (38) is not valid in that casef]. 



TYPE IIB SUPERSTRING FROM THE NBI MATRIX MODEL 

We can now summarize our results in the following rather long formula, 

Z = J dA^difjdY exp f^Tr(y -1 [A M ,A,] 2 ) - (3 TrY - (n - i) TrlnF 

= const. J dA^dipdM ^det yj '-[A^, ~Ajf\ 

xexp (-Tr ((^3 + M 2 ) A v \ A - |Tr (^[A^])^ 

-> const. J dX^dipdM ^det ^{X^X u } 2 ^j 

xexp(- J dadr ((y/rf + M 2 )y / {X^Z} 2 - y ^})J J , (39) 



where the last expression is valid for n —>■ oo. In this expression the fields ip and M 
have expansions similar to eqs. (21) and (32), and some normalization constants have 
been absorbed in a and (3 in the last formula above. 

The functional integration over M in eq. (39) is just Gaussian and can of course 
easily be performed, 



const. / dX^dip 



det y/{X M ,Xj 2 X V2 
Vety/{X^X u y 

xexpf-J dadr \J^^{X^X V ¥ - j(^{X^})\\ . (40) 

The two determinants in this expression arises from different types of Gaussian inte- 
grations, the "det" beeing defined through eq. (18) and the subsequent limit n — > oo, 
whereas the u T>et n determinant comes from the continuum integral over M. Naively 
one would tend to identify these two determinants, so that the fraction containing them 
is just one, 

det y/{X„X^ ^ i 

If so, the NBI action gives exactly the Nambu-Goto version of the Green-Schwarz type 
IIB superstring. 

However, Zarembo 9 has pointed out to me that the situation can be more com- 
plicated. For example, from eq. (18) one sees that the "det" determinant (= detz for 



^For such strings where the end points are actually separated by a 7arge distance, the interchange of 
limits in (34) is probably allowed. Thus at large distances the string picture is most likely right for 
the NBI matrix model even without supersymmetry. At shorter distances near the critical one, this 
picture breaks down, and the sine function in (31) cannot be approximated by its first term in a 
power series expansion. 
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n — > oo) is subdominant relative to the factor Yl(zi + Zk) occurring in eq. (18), and 
hence can be ignored in the limit n — > oo. In this case, one has instead of (41) 

If so, the Vet determinant survives in the measure. However, even if this is so, this factor 
is rather harmless: correlation functions are invariant, since the measure is multiplied 
by a constant factor under reparametrizations. This factor does not depend on the 
fields and cancels in the correlation functions 9 . 

Perhaps the right answer depends on how exactly the continuum limit is con- 
structed, because in order to interpret Vet a regulator is needed. 

It should also be mentioned that Chekhov and Zarembo 10 have discussed models 
somewhat different from the NBI model, and have also discussed the measure in more 
details. 



A SADDLE POINT AND THE VIRTUAL EULER NUMBER 

We shall now study the saddle point of the NBI action. By variation of the 
A^— fields we obtain the classical equation of motion 

[A ll ,{Y- 1 ,[A lt ,A v ]}]=0. (43) 
This equation was studied by Kristjansen and me 11 . The solution is 

[A„A v ]*=imr,Y}, (44) 

where m^ u is a matrix with repect to the space indices. In the saddle point the action 
has the value 

^Nm dle = (P + m%,a/4) TrF + (n - 1/2) Tr In Y. (45) 

In order to have a non-trivial n — > oo limit, it is necessary that a and (5 are of order n. 
It should be stressed that this does not imply the usual classical limit in string theory, 
as explained in details in ref. 11. 

In addition to the terms exhibited above, there are of course subdominant terms 
arising from the expansion of around the classical solution. These terms are ignored 
in the following. Therefore, at the A M -saddle point we have the integral 11 (a/n and 
P/n are both of order one) 

Z saddi c = J dYexp[-n{(P + ml u a/4:)/nTrY + t Trlny}], with t = 1 - l/2n. (46) 

This functional integral is of the Penner type 12 . For the value of the parameter t needed 
in the saddle point, the F-integral actually diverges. However, by analytic continuation 
a la the gamma function for negative argument one can start by defining the integral 
for negative t, and then ultimately continue back to positive t. It the turns out that 
t — 1 is a critical point, and in the vicinity of this point one can define a double scaling 
limit with the "cosmological constant" 

H = (1 - t)n = fixed. (47) 

We see that with the value of t given by eq. (46), /x =l/2=fixed quite automatically! 
Thus, we do not need to make any special assumptions in order to have this double 
scaling limit in the NBI model. 
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What is the meaning of the Penner model in the double scaling limit? An asymp- 
totic expansion in fi can be made. Consider the "free energy" F, Z sa,ddle = e F , then 

oo 

F(fi) = F (fi) + F 1 ( f i) + Y,xy- 29 , ^ = 1/2- (48) 

9=2 

Here Xg is the "virtual Euler number" for moduli space of Riemann surfaces with genus 
g, which is well known to be relevant for strings. One has 12 

Xa = 20(20-2)' (49) 

where B 2g are the Bernoulli numbers. These have positive sign, and blow up factorially, 
so the sum defining F is not Borel summable. This is also well known to be the case 
for genus expansion of string theories. 

The physical interpretation of this result is that the field Y captures the Euler 
characteristic of moduli space of Riemann surfaces. Therefore it is quite likely that 
the NBI model encodes non-perturbative information on Riemann surfaces generated 
by moduli space. It should be remembered that in string theories one usually sum the 
functional integral over g, however here this seems to be already included. It must 
be admitted that the virtual Euler number represents very global properties of moduli 
space, and certainly more details are needed before one can claim a good understanding 
of the non-perturbative nature of this model. 

Recently Soloviev 13 has commented on "a curious relation" between Siegel's model 14 
of random lattice strings and the above saddle point approximation to the NBI model. 
This comes about if one starts from Siegel's T-self-dual matrix model 

S = Tr Q$ 2 + n ln(l - gty^J , (50) 

where $ is a Hermitean n x n matrix and g is a constant. It was then pointed out by 
Soloviev 13 that if one makes the substitution 

gY=l-g$, (51) 

and perform the limit n — » oo, g — > 0, gn =fixed, then one obtains 

S — > n Tr (const. Y + n In Y) + irrelevant const.. (52) 

This is, however, precisely the saddle point expression (46) for the NBI model. This 
sadle point is therefore a weak string coupling limit of the Siegel matrix model 13 . For 
arbitrary coupling there is, however, an additional F 2 -term in the Siegel action, and 
hence it was suggested that perhaps the potential (3) should have an additional TrY 2 
term 13 . Of course, a similar statement can be made about the NBI model, where there 
are various corrections to the saddle point expansion. 
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